Electronic band gap and transport in Fibonacci quasi-periodic graphene superlattice 
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We investigate electronic band gap and transport in Fibonacci quasi-periodic graphene superlat- 
tice. It is found that such structure can possess a zero-fc gap which exists in all Fibonacci sequences. 
Different from Bragg gap, zero-fc gap associated with Dirac point is less sensitive to the incidence 
angle and lattice constants. The defect mode appeared inside the zero-fc gap has a great effect on 
transmission, conductance and shot noise, which can be applicable to control the electron transport. 



Graphene, a monolayer of carbon atoms tightly packed 
into a honeycomb lattice, has attracted great interest 
in graphene-based nanoelectronic and optoelectronic de- 
vices [l||, since it was fabricated by Novoselov and Geim 
et al. in 2004 ■ In graphene, the unique band structure 
with the valance and conduction bands touching at Dirac 
point (DP) leads to the fact that electrons around the 
Fermi level can be described as the masslcss rclativistic 
Dirac electrons, resulting in the linear energy dispersion 
relation. As a consequence, there are a great number 
of electronic properties, such as the half-integer quan- 
tum Hall effect , the minimum conductivity Q , and 
Klein tunneling [6j. In particular, Klein tunneling and 
perfect transmission are crucial for electron transport in 
various graphene heterostructures 0], i.e. single barrier 
H and n-p-n junctions Q . 

Motivated by the ex per imental realization of graphene 
superlattice (GSL) [lX>jl2l ]. electronic bandgap structures 
and transport properties in GSLs with electrostatic po- 
tential and magnetic barrier have been extensively in- 
vestigated [l3|-|22( , since the conventional semiconductor 
superlattices are successful in controlling the electronic 
structures and the extension to graphene may give rise to 
different features and applications. For instance, DP ap- 
pears in the GSL [l4|, [H[ , and it is exactly located at the 
energy with the zero-fc gap pTj . Interestingly, the zero-fc 
gap associated with DP is insensitive to the lattice pa- 
rameter changes in contrast with the behavior exhibited 
by Bragg gaps fl7j . This gap is analogous to photonic 
zero-n gap in the photonic crystals containing negative- 
index and positive- index materials [20| . and originates 
from a zero total phase [23| . Accordingly, the zero-fc gap 
is robust against the lattice constants, structural disorder 
[l7j|, and external magnetic field 18 1, and thus is better 
to control the electron transport in GSL. 

In this Letter, we will investigate electronic band gap 
and transport in Fibonacci quasi-periodic GSLs in the 
fashion analogous to photonic crystal with metamatcrials 
j23l - |25j . As we know, the quasi-periodic GSL is classified 
as intermediate between ordered and disordered systems 
fl9l [20| , which has significant and common features like 
fractal spectrum and self-similar behavior [2l|, [22[ • How- 
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FIG. 1. (Color online) (a) Example of quasi-periodic mono- 
layer GSL, corresponding to Fibonacci sequence S3, (b) The 
schematic profiles of the potentials Va and Vb- 



Corresponding author. Email: xchcn@shu.cdu.cn 



ever, what we concentrate on here is the electronic band 
gap and DP in such quasi-periodic system. We find that 
zero-fc gap happens in all Fibonacci sequences, which re- 
sults in the robust transmission properties, conductance 
and shot noise at the DP. 

Consider quasi-periodic monolayer GSLs with the 
structure in each cell following the Fibonacci sequence, 
Sj, by a recurrent relation Sj+i = {Sj, Sj-i}, with 
Sq = {B} and Si = {A} with j is the generation num- 
ber of the Fibonacci unit cell, the first few sequences 
are S* 2 = {AB}, S 3 = {ABA}, S 4 = {ABAAB} and 
so on. Elements A and B are considered as the alter- 
nating barriers Va and wells Vb with the width d,A and 
ds, respectively. As an example, the third-generation 
Fibonacci structure (ABA)" 1 with the number of peri- 
ods, m, is shown in Fig. [TJ Generally, in the vicinity 
of the K point and in the presence of a potential V(x), 
the charge carriers are described by the Dirac-like equa- 
tion, H = —ihvFO' ■ V + V(x) where the Fermi velocity 
vp w 10 6 m/s, and a = (a x ,a y ) are the Pauli matri- 
ces. Due to the translation invariance in the y direc- 
tion, the solution of above equation for a given incident 
energy E and potential barrier Vj can be presented as 



FIG. 2. (Color online) Energy band (a) and transmission 
spectrum (b) for the Fibonacci quasi-periodic GSLs with 
d B /d A = 1, where V A = 50 meV, V B = meV, 6 = 20°, 
m = 16, (a): d,A = 20 nm; (b): d A = 15 nm (dotted black 
line), <1a — 20 nm (solid red line), and cIa = 25 nm (dashed 
blue line). 
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amplitude of the forward (backward) propagating wave 
The wave functions at any two positions x and x + Ax 
inside the jth potential can be related via the transfer 
matrix 171: 
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with Qj = arcsin(fcj,/fcj) 
coefheient t 



As a result, the transmission 
t(E, k y ) is found to be 



2 cos 9q 



m 2 i 



(2) 



where 8 and t are incidence and exit angles (see Fig. 
[1]), and m,ij(i,j = 1,2) is the matrix element of total 
transfer matrix, Xn = n^Li ^-h connecting the incident 
and exit ends, and N is the total number of layers of the 
grapheme supcrlattice. Once the transmission coefficient 
is obtained, the total conductance G of the system at zero 
temperature is given as follows, G = G J^ 2 T cos 9 d9 
[HI, where T = \t\ 2 and G = 2e 2 mv F L y /h 2 and L y 
is the width of the graphene stripe in the y direction. 
Meanwhile, the Fano factor is given by F = jZ„/2 T(l — 

T) cos 9 d6 1 j^ 2 /2 T cos o d0 o [13] . 

Fig. [2] shows the energy bands and transmission spec- 
trum in various Fibonacci quasi-periodic GSLs (i.e., from 
5-2 t° St). Besides the distribution of energy bands like 
Cantor- like set j22| , what we have discovered here is that 
the zero-fc gaps exist in all Fibonacci levels. In Fig. [2] (a), 
there are several broad forbidden gaps opened for each 
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FIG. 3. (Color online) Transmission spectrum for the Fi- 
bonacci quasi-periodic GSL, (ABA) m , where m = 16, (a) 
d,A = 20 nm, (b) d,A = 25 nm, and the other parameters are 
the same as those in Fig. [21 (c) and (d) are the electronic 
band structures corresponding to (a) and (b). 



Fibonacci level in the considered energy range. Among 
these forbidden gaps, we notice that the position and size 
of zero-fc gaps are almost robust against the Fibonacci 
levels. In fact, the Fibonacci structure S2 is exactly 
the GSL, (AB) m . The condition for zero-fc gap is given 
by qA&A = — qB^B at 6a = 0, which provides the DP, 
E = Va/(1 + (1b /(1a), for the special case of Va 7^ and 
Vb =_0 O- For the higher Fibonacci level S3 to 5*7, the 
zero-fc gaps become stabilized with the fixed position and 
size, although the location of zero-k gaps is slightly dif- 
ferent from that for Fibonacci level £2. Furthermore, we 
demonstrate, in Fig. H](b), that such gap depends only 
on the ratio of lattice constants, and is insensitive to the 
lattice parameters. On the contrary, the position and size 
of Bragg gaps in a higher energy range change sensitively 
with the Fibonacci level and lattice parameters. 

Fig. [3] (a) and (b) further display the influences of 
the incidence angle on the transmission spectrum in the 
Fibonacci quasi-periodic GSL, (AB) m , corresponding to 
Fibonacci sequence S3. It is apparent that the zero-k 
gap is independent of the lattice constants and is weakly 
dependent on the incidence angle. To understand better, 
the electronic dispersion at any incidence angle, based on 
the Bloch's theorem, is written as, 



= cos (2q A dA) cos {qBd B ) 
sin (2q A dA) sin {q B d B ), (3) 



cos (&.A) = -Tv[M A M B M A ] 

sin 6 a sin 6b — 1 
cos 6a cos 6b 

where A = 2dA + g?b is the length of the unit cell. 
Therefore, the location of the DP is given by 2^6?^ = 
~qs (1b at 6a — 0. For the structure considered here, 
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FIG. 4. (Color online) Conductance (a) and Fano factor (b) as 
a function of Fermi energy in Fibonacci quasi-periodic GSL, 
(ABA) 1 ", where m = 16, d A = 20 nm (dotted blue line) and 
ds = 25 nm (solid red line), and the other parameters are the 
same as those in Fig. [2] 
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FIG. 5. (Color online) Effect of defect mode on the trans- 
mission (a), conductance (b), Fano factor (c) in GSL, where 
the defect do = 80 nm, Vn = 55 meV, d A = 20 nm, and the 
other parameters are the same as those in Fig. [2] 



(Va 7^ and Vg = 0), the DP is exactly located at 
E = Va/(1 + (Ib/ScIa), which means that the zero-fc gap 
depends only on the ratio, cLb/cLa, instead of d,A and ds 
themselves. Fig. [3] (c) and (d) show that a band gap 
opens at E = 33.3 meV, which is different from E = 25 
meV for Fibonacci sequence 5*2 , in which ds I Aa = 1 and 
Va = 50 meV. In fact, DP for other Fibonacci sequences 
can be further calculated as E = Va/W + ds/TjdA], 
where Tj is the ratio of numbers of layer A and B, and 
lim^oo Tj = (1 + V5)/2 = 1.618. 

Fig. 2] shows the conductance and the Fano fac- 
tor with the difference lattice constants. Remarkably, 
the angular-averaged conductance reaches the minimum 
value at the DP, while the Fano factor exists a peak in 
the vicinity of DP with value approximately F = 1/3 
[l8l . [27j . The conductance and the Fano factor shows the 
robust properties, since the DP does not shift, when the 
lattice constants dA and ds are changed simultaneously. 
In Fig. [5J we further shed light on the effect of localized 



defect mode in GSL, (ABA) 8 D{ABA) 8 , where the defect 
layer D with do = 80 nm and Vb = 55 meV. Compared 
to that inside the Bragg gap, the defect mode inside the 
zero-fc gap remains almost invariant with the incidence 
angles and lattice constants. Due to the existence of the 
defect mode, the conductance is greatly enhanced, while 
the Fano factor is strongly suppressed, as shown in Fig. 
[5] This suggests that the electron transport can be mod- 
ulated by the defect mode. 

In summary, using the transfer matrix method, we 
have investigated the electronic band gap and transport 
in the Fibonacci quasi-periodic GSL. It is shown that 
the zero-fc gap and the defect mode are robust against 
the lattice constants and incidence angle, which is useful 
to control electron transport. We hope such Fibonacci 
structure will have applications in graphene-based elec- 
tronic omnidirectional reflector and filters. 

This work was supported by the NSFC (Grant Nos. 
60806041 and 61176118), and the Shanghai Leading Aca- 
demic Discipline Program (Grant No. S30105). X. C. 
also acknowledges Juan de la Cierva Programme, the 
Basque Government (Grant No. IT472-10) and MICINN 
(Grant No. FIS2009-12773-C02-01). 



[1] A. H. Castro Neto, F. Guinea, N. M. R. Peres, K. S. 

Novoselov, A. K. Geim, Rev. Mod. Phys. 81, 109 (2009). 
[2] K. S. Novoselov, A. K. Geim, S. V. Morozov, D. Jiang, 

Y. Zhang, S. V. Dubonos, I. V. Grigorieva, A. A. Firsov, 

Science 306, 666 (2004). 
[3] K. S. Novoselov, A. K. Geim, S. V. Morozov, D. Jiang, 

M. I. Katsnelson, I. V. Grigorieva, S. V. Dubonons, and 

A. A. Firsov, Nature (London) 438, 197 (2005). 
[4] Y. Zhang Y. W. Tan, H. L. Stromer, and P. Kim, Nature 

(London) 438, 201 (2005). 
[5] V. P. Gusynin, and S. G. Shararpov, Phys. Rev. Lett. 

95, 146801 (2005). 
[6] M. I. Katsnelson, K. S. Novoselov, and A. K. Geim, Nat. 

Phys. 2, 620 (2006). 
[7] A. F. Young and P. Kim, Annu. Rev. Condens. Matter 



Phys. 2, 101 (2011). 
[8] X. Chen and J.-W. Tao, Appl. Phys. Lett. 94, 262102 
(2009). 

[9] V. V. Cheianov and V. I. Fal'ko, Phys. Rev. B 74, 

041403(R) (2006). 
[10] J. C. Meyer, C. O. Girit, M. F. Crommie, and A. Zettl, 

Appl. Phys. Lett. 92, 123110 (2008). 
[11] S. Marchini, S. Giinther, and J. Wintterlin, Phys. Rev. 

B 76, 075429 (2007). 
[12] A. L. Vazquez de Parga, F. Calleja, B. Borca, M. C. G. 

Passeggi, Jr., J. J. Hinarejos, F. Guinea, and R. Miranda, 

Phys. Rev. Lett. 100, 056807 (2008). 
[13] C.-X. Bai and X.-D. Zhang, Phys. Rev. B 76, 075430 

(2007). 

[14] M. Barbier, F. M. Peeters, and P. Vasilopoulos, Phys. 



4 



Rev. B 80, 205415 (2009); 81, 075438 (2010). 
[15] L. Brey and H. A. Fertig, Phys. Rev. Lett. 103, 046809 

(2009) . 

[16] C. H. Park, L. Yang, Y. W. Son, M. L. Cohen, and S. G. 

Louie, Phys. Rev. Lett. 101, 126804 (2008). 
[17] L.-G. Wang and S.-Y. Zhu, Phys. Rev. B 81, 205444 

(2010) ; L.-G. Wang and X. Chen, J. Appl. Phys. 109, 
033710 (2011). 

[18] X.-X. Guo, D. Liu, and Y.-X. Li, Appl. Phys. Lett. 98, 
242101 (2011). 

[19] N. Abedpour, A. Esmailpour, R. Asgari, and M. R. R. 
Tabar, Phys. Rev. B 79, 165412 (2009). 

[20] Y. P. Bliokh, V. Freilikher, S. Savel'ev, and F. Nori, Phys. 
Rev. B 79, 075123 (2009). 

[21] S. Mukhopadhyay, R. Biswas, and C. Sinha, Phys. Sta- 
tus. Solidi (b) 247, 342 (2009). 



[22] S. H. R. Sena, J. M. Pereira Jr, G. A. Farias, M. S. 

Vasconcelos, and E. L. Albuquerque, J. Phys.: Condens. 

Matter. 22, 465305 (2010). 
[23] J. Li, L. Zhou, C. T. Chan, and P. Sheng, Phys. Rev. 

Lett. 90, 083901 (2003). 
[24] A. Bruno- Alfonoso, E. Reyes-Gomez, S. B. Cavalcanti, 

and L. E. Oliveira, Phys. Rev. A 78, 035801 (2008). 
[25] L.-W. Zhang, K. Fang, G.-Q. Du, H.-T. Jiang, and J.-F. 

Zhao, Opt. Commun. 284, 703 (2011). 
[26] S. Datta, Electronic Transport in Mesoscopic Sys- 
tems (Cambridge University Press, Cambridge, England, 

1995). 

[27] J. Tworzydlo, B. Trauzettel. M. Titov, A. Rycerz, and 
C. W. J. Beenakker, Phys. Rev. Lett. 96, 246802 (2006). 



